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ABSTRACT 



Aims. We present an analysis of the Lyman a forests of five quasar spectra in the near UV. Properties of the intergalactic medium 
(IGM) at an intermediate redshift interval (0.9 < z < 1.9) are studied. The amount of baryons in the diffuse photoionised IGM 
and the warm-hot intergalactic medium (WHIM) are traced to get constraints on the redshift evolution of the different phases of the 
intergalactic gas. 

Methods. The baryon density of the diffuse IGM is determined with photoionisation calculations under the assumption of local 
hydrostatic equilibrium. We assume that the gas is ionised by a metagalactic background radiation with a Haardt & Madau (2001) 
spectrum. The WHIM is traced with broad Lyman a (BLA) absorption. The properties of a number of BLA detections are studied. 
Under the assumption of collisional ionisation equilibrium a lower limit to the baryon density could be estimated. 
Results. It is found that the diffuse photoionised IGM contains at least ~ 25% of the total baryonic matter at redshifts 1 < z < 2. 
For the WHIM a lower limit of ~ 2.4% could be determined. Furthermore the data indicates that the intergalactic gas is in a state of 
evolution at z ~ 1.5. We confirm that a considerable part of the WHIM is created between z = 1 and z = 2. 

Key words, cosmology: observations - quasars: absorption lines - intergalactic medium 



1. Introduction 

The baryonic matter in the universe can be divided into four 
gas phases. The largest fraction of baryons resides in the dif- 
fuse IGM which is photoionised by a metagalactic UV back- 
ground radiation. While at higher redshifts z > 3 about 90% 
of the baryonic matte r is assembled in the photoionised IGM 
dWeinberget al.lll997h . in the local universe this g as phase in- 
herits a much smaller fraction of matter (~ 30%) (Penton et al.' 
2004; Lehner et al. 2007). This incident is called the missing 
baryon problem. It is now widely believed that the gas has been 
shock-heated by gravitational collapse to higher temperature s 
r ~ 10^ - 10^ K thus forming the WHIM C Valageas et alll2002l) . 
It is expected to contribute a major fraction of the baryonic 
matter at lo w redshifts (Cen & Ostriker 1999; Dave et al. 2001; 
lDavel2063h . Hot gas that is gravitationally bound in galaxy clus- 
ters as well as condensed objects like stars and cool galactic gas 
only contain a bout 10% of the total baryons each (Dave 2003; 
iBell et al.ll2003h . 

The WHIM is a low density («h ~ 10"*' - lO^^^cm^^) high 
temperature plasma and is thus hard to detect. Current X-Ray 
observatories are not sensitive enough to detect its diffuse emis- 
sion which is suspected to have a very low surface brightness. 
However, the analysis of quasar absorption features seems to be 
a promising approach. 

Oxygen is a favoured element since it is relatively abundant 
and its transitions have large oscillator strengths. The O vi ion 
traces gas at a temperature of T ~ 3 ■ 10^ K, characteristic for 
the WHIM. At low redshifts (z < 0.5) various detections of 
Ovi absorbers that are associated with the WHIM have been 



iDanforth & Shulll 120051) . These measurements imply a baryon 
density of the WHIM in Ovi absorption of Qb(Ovi) > 
0.0022/!7(j, assuming that 20 percent of the oxygen is fivefold 
ionised and that the mean oxygen abundance of the IGM is 0. 1 
solar. This relatively low value is explained with the absence of 
O VI in hotter gas with temperatures T > 10^ K which might con- 
tain a large fraction of baryons. Higher oxygen transition (O vii 
and O viii) could be detected with X-ray absorption measure- 
ments, however, the low spectral resolution of cuiTent X-ray ob- 
servatories limits the possibilities in this research field. 



Currently the search for broad Lyman a (BLA) absorption 
of neutral hydrogen seems to be the most promising approach. 
Due to the high gas temperature, WHIM absorption lines will 
be strongly broadened and due to the expected small neutral 
fraction (/hi < 10"^) most lines will be shallow. Some de- 
tections o f BLAs at low redshifts have been reported in the 
last years ("Tripp et al."200r; Richter et al. "2004"; Sembach et al 
2004; Richter etal. 2006b; Lehner etal. 2007). Richter etal 
(2006b) derived a lower limit of Ob(BLA) = 0.0027h-i for the 



reported (Tripp et al. 2000; Oegerle et al. 2000; Tripp & Sayag 
2000; Chen & Prochaska 2000;' Savage et al. 2002; Richter et al. 
2004 rSembach et aL ,2004 .Savage et al.i ,2005.: ,Tripp ,2005c 



WHIM baryon density. It is expected that the WH IM makes up 
at lea st 20% of the baryons in the local universe dLehner et al.l 
12001 . 



In this paper we estimate the amount of baryons residing in 
the IGM, and the WHIM in particular, at an intermediate red- 
shift interval (0.9 < z < 1.9). The evolution of the gas phases 
is of particular interest. In Sect. |2] some properties of the IGM 
are discussed. Sect.[3]describes the data used in this work. BLA 
detections in the sightline of five quasars are reported and dis- 
cussed in Sect.m A summary of the results will be given in Sect. 
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2. Properties of the IGM 

2.1. The baryon density 

To make an estimation of the baryon content of the IGM, the 
total hydrogen density has to be known. Since the analysis of 
quasar absorption lines just gives the column density of neutral 
hydrogen and the gas is highly ionised, the ionisation fraction is 
an important factor for an estimation of the baryon density. With 
knowledge of the total hydrogen column density A^h of each ab- 
sorber, the density parameter Q can be summed up: 



PcC 



(1) 



where ^ - 1 .3 is a correction factor for the helium abundance, 
niH the hydrogen mass, Hq the Hubble constant, pc the critical 
density, c the speed of light and AX is the comoving path length 
(Eq.|2]|. The index / denotes a Lyman a system along the line of 
sight 

The comoving path length AX was calculated using a ACDM 
cosmology: 



dX- 



l+z 



VnM(l+z) + QA/(l+z) 



dz. 



(2) 



with the dark energy contribution Qa = 0.73, the dark matter 
density Qqm — 0.23 and the baryonic matter density Qbm = 0.04 
jSpergel et al.ll2003[) . 

Since the diffuse photoionised IGM and the WHIM are gov- 
erned by different physics, different approaches have to be made 
to determine the ionisation fraction. 

2.2. Photoionised IGM 

The ionisation fraction of the diffuse photoionised IGM is 
mainly dependent on the photoionisation rate of the ionising 
background radiation, the gas temperature and the gas density. 
In this section some properties of the photoionised IGM are dis- 
cussed that will be used to estimate the density parameter. 

The neutral fraction of a highly ionised optically thin gas can 
be written as 



»HI 



(3) 



where q-hh ~ 4 ■ lO^'^'T""-^^ cm-'s"' is the recombination co- 
efficient of ionised hydrogen, T is the gas temperature and 
r - ■ lO'^s ' is the photoionisation rate of the UV back- 
ground radiation. 

Under the assumption of local hydrostatic equilibrium, the 
typical absorber size should be of the same order as the local 
Jeans length. In this case the total hydrogen column density A^hj 
can be expressed in terms of the total hydr ogen volume d ensity 
riii and the gas temperature T - T4 ■ 10* K ( ISchavell200ll) : 



H,J 



1.6 



in21 -2 1/2,7,1/2 

10 cm T^' 



A. 

0.16 



(4) 



where fg is the fraction of mass in gas. On the scales of interest fg 
will not be far from its universal value f„ - ~ 0.16, however 
in cold collapsed clumps it could be close to unity. 

Substituting the neutral fraction (Eq. |3]l into Eq.|4] gives the 
neutral hydrogen column density A^hi: 



A^Hi ~ 2.3 ■ 10'^ cm"^ 



1 10-5 cm-^l 



7,-0.261—1 1 fs 



(5) 



The photoionisation rate F can be integrated with knowledge 
of the UV background intensity J{v) and the photoionisation 
cross section a(v): 



r 



= 4. f:^., 



(v)dv. 



(6) 



where h is the Planck constant and vq is the frequency at the 
hydrogen ionisation energy. 

2.3. WHIM 

The WHIM is a low density, high temperature plasma and is 
thus expected to create broad and shallow absorption features in 
neutral hydrogen . The ionisation fraction ca n be estimated, fol- 
lowing a work bv Sutherland & Dopital (Il993l) . They have inves- 
tigated cooling functions for low-density astrophysical plasmas 
under the assumption of collisional ionisation equilibrium. For 
the temperature range of interest (T ~ 10^ - 10^ K) the ionisa- 



tion fraction fio 



polynomial (iRichter et aTjl2006bh : " 
log/ionCD ^ -13.9 + 5.4 log r - 0.33(logr)2 



Tp can be approximated as a 



(7) 



It is solely dependent on the gas temperature T and does not 
vary over a wide range of metallicities. Howev er, it is just valid 
for pure collisional ionisation equilibrium. Ric hter et alJ (l2006ah 
found that photoionisation is important for the WHIM as well 
and that neglecting it might underestimate the baryon density up 
to 50%. The best fit to their simulated data is a linear relation 
between the ionisation fraction and the gas temperature: 



log/ion(r)* -0.75 + 1.25 log r. 



(8) 



The temperature of an absorber can be estimated from the 
measured Doppler parameter. The broadening of the line shape 
is created by a combination of thermal and non-thermal effects: 

b = -^^^i'^^'^^on-th' ^'^^ P™^ thermal broadening the Doppler pa- 
rameter b would be directly connected with the gas temperature 
T: 



b,h - 



2kT 
m 



(9) 



where k is the Boltzmann constant and m is the particle mass. 
WHIM H I absorbers with temperatures T > 10^ K will thus cre- 
ate absorption features with Doppler parameters b > 40 km s"'. 

A determination of the gas temperature requires a de- 
tail ed study of the l ine shape of each absorber. According to 
Ric hter et al.1 (l2006bh non-thermal broadening effects should be 
detectable in high resolution spectra with a good signal to noise 
ratio (S /N). The most important non-thermal broadening mech- 
anisms are macroscopic velocities and turbulence within the gas. 
These should, in general, produce non-Gaussian line shapes. 
Though there could be cases in which macroscopic velocities 
create a Gaussian shaped broadening, e.g. a symmetric infall 
of gas towards the line of sight, these should not be dominant. 
Since WHIM absorbers will have Doppler parameters of at least 
b > 40kms-' the Hubble flow, which is suspected to contribute 
with ~ 7 km s ' per hundred kpc absorber width, should not be 
dominant for typical absorber sizes. 

Note that the ionisation fraction increases with the gas 
temperature. Consequently absorbers with higher temperatures 
would inhabit more baryons than colder absorbers with the 
same column density. Therefore systems that inherit the high- 
est amount of baryons are the most broadened and might thus be 
the hardest to detect. 
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Table 1. Quasar overview 



QSO 


z 


SIN 


A.1 [A] 


Az 


Azbl^ 


HE 0515-4414 


1.73 


10-50 


2310-3260 


0.782 


0.240 


HE 0141-3932 


1.80 


-25 


3061-3384 


0.266 


0.101 


HE 2225-2258 


1.89 


-25 


3057-3478 


0.346 


0.033 


HE 0429-4901 


1.94 


- 8 


3188-3538 


0.288 


0.126 


PG 1634+706 


1.34 


5-50 


2310-2790 


0.395 


0.255 



" Just the regions with {S IN) > 15 that are not too much populated 
with Lyman a forest lines could be used. 



3. Observations 

The Lyman a forests of five quasars were analysed in the near 
UV. A wavelength range from 2300 A to 3500 A was covered. 
Up to the atmospheric cut-off at ~ 3000 A the data was taken 
with the Space Telescope Imaging Spectrograph (STIS) on board 
the Hubble Space Telescope. Higher wavelengths were covered 
with the UV Echelle Spectrograph ( UVES) at the ESQ Very 
Large Telescope. For HE 0515-4414 jReimers et al.ll 19981) data 
were available from both instruments thus covering a wide wave- 
length range. The E230M grating of STIS gives a resolution of 
R ~ 30000. The UVES data were taken with a slit width of 
0.8" resulting in a resolution of R ^ 50000. The spectra of 
three additional quasars from the Hamburg ESO survey were 
used, that have been obtained in the ESO service mode for pre- 
vious studies. HE 0141-3932 and HE 2225-2258 were observed 
with a slit width of 1 .0" resulting in a resolution of R ^ 40000 
dWisotzki et al.ll2000l) and HE 0429-4901 again with a slit width 
of 0.8". Furthermore da ta from the STIS archive was used. The 
data for PG 1634H-706 jSchmidt & Greenlll983h were taken in 
two programs (S. Buries 7292/7 and B. Januzi 8312/8) also us- 
ing the E230M grating. 

Exce pt for HE 0515-4414, pipeline data was used for all 
quasars. iJanknecht et al.l(l2006bh have reworked the reduction of 
HE 0515-4414 and have done vacuum and barycentric coiTec- 
tions to the wavelength scale as well as a co-addition of multiple 
exposures for all quasars. Table [3] shows the redshift z, the sig- 
nal to noise ratio S IN per pixel, the used wavelength range IS.A 
with the corresponding redshift range Az and the effective red- 
shift range for the search for BLAs Azbla- 

A complete list o f all Lyman a forest absorbers was taken 
from IJanknecht et al.l (|2006a). It was used for the photoionisa- 
tion calculations and as a basis for the search for BLA absorbers. 
The S IN varies strongly over the spectrum of each quasar. Since 
the possibility of detecting broad and shallow features is strongly 
dependent on the noise, not the whole available wavelength 
range was usable for the search for BLAs. An average signal 
to noise ratio of {S IN) > 15 was used as a lower limit to the 
data quality. Furthermore a region of Av = 5 000 km s"' around 
the Lyman a emission line of each quasar was excluded to avoid 
influences of the proximity effect. This resulted in a usable red- 
shift range of Azbla = 0.723. The STIS and the UVES data 
cover nearly the same range of analysed data. 

4. Results 

4.1. Photoionisation Calculations 

We used lHaardt & Madaul ( 1200 lb spectra (HM) to calculate the 
photoionisation rates F at different redshifts (Eq. |6]l. For the 
mean intensity of the ionising background at the hydrogen ion- 



isation edge the values provided by HM were used. They in- 
crease from log Jy = -21.61 for z ~ 0.9 to log Jy = -21.29 for 
z ~ 1.9. The fraction of mass in gas was set to its universal value 
/g-0.16. 

The gas temperature is supposed to be about T ~ 
lO'^K, however, a more precise determination is very difficult. 
Therefore, in the following considerations the gas temperature 
will be left variable. Using Eq. |5] the neutral hydrogen volume 
density can be approximated for each absorber. The resulting to- 
tal hydrogen column densities (Eq. |4]l can be summe d up to a 
total baryon density using Eq.\t\ Since ISchavd (1200 Ih restricts 
his analysis to optically thin clouds, only absorbers with column 
densities logA^Hi < 17 were used. A relatively low significance 
level of SL = ^ > 1, where W is the observed equivalent 
width and crw its Icr error, was used for the selection of Lyman 
a lines ( Janknecht et al. 2006b). The resulting sample contains 
624 Lyman a absorbers in the lines of sight of five quasars, cov- 
ering a total redshift interval of Az - 2.08, corresponding to a 
comoving path length of AX - 5.465. 

The complete sample results in a baryon density of 
Qb(LyQ') = 0.010 + 0.001 ■ T^^''. Dividing the sample into two 
redshift intervals shows an observable increase of the baryon 
content with redshift. While for 0.9 < z < 1.5 the density 
parameter is found to be n^(Lya) = 0.008 ± 0.001 ■ 
at higher redshifts 1.5 > z > 1.9 it increases to D.b{Lya) = 
0.012 ± 0.001 ■ r^'^^. The errors are propagated from the er- 
rors of the column density fits. The real uncertainty would be 
much higher. It should be remembered that the absolute value 
of the baryon density depends directly on the mean intensity Jy 
of the ionising background, which is not well known. A change 
in the photoionisation rate F of 20% coul d cause a change i n 
the derived baryon density of about 10% (iLehner et al.ll2007l) . 
However, the relative numbers which trace a possible variation 
of the baryon density would not be sensitive to this error source. 
Note that the results here are just applicable for a constant gas 
temperature. Since the temperature of the ph otoionised gas is ex- 
pected to decrease with decreasing redshift dRicotti et al.ll200()t 
'Schave et al."2000) the evolution could be even stronger 

Schave (2001) proposes a more statistical approach to the 
problem. In the literature observational results are often ex- 
pressed in terms of the number of absorption lines per unit ab- 
sorption distance X and per unit column density /(A^hi.z)- Using 
this function, the density parameter can be integrated: 

Qb = ^^^^ NH,—fiNnuZ)dNn,. (10) 

PcC Jn,„,„ "Hi 

Combining this equation with equations[3]and|5]yields 
Qb~2.2- lO-V^Fj^^rf^- Nll^f{NH„z)dNH„ (11) 

where h = lookm?*! Mpc-' ■ '^he fraction of mas s in gas fg was 
again s et to its universal value fg = 0.16. IJanknecht et al.l 
(l2006bl) found that the column density distribution can be ex- 
pressed as a single power law, /(Nhi) ~ AN^f^ with /3 = 
1.60 ±0.03 and log A = 9.4 + 0.3 over the total observed redshift 
range. Adopting an averaged photoionisation rate F12 - 1.07, 
according to z ~ 1.5, and the integration limits Nmin = 12.9 and 
A^max = 15.7 dJanknecht et a l. 2006a) results in a baryon den- 
sity of Qb(LyQ;) = 0.009 ■ 7'^'^'' which is within the Icr error 
interval of the previous result. However, the density parameter 
strongl y depends on the slope of the column density distribu- 
tion jS. IJanknecht et alj (l2006bl) found a slightly different slope 



4 



N. Prause et al.: The baryon density at z = 0.9 - 1.9 



for the higher redshift sample 1.5 < z < 2 of yS = 1.55 ± 0.04. 
For this analysis the uncertainty in the determination of (3 has 
a very big impact. Using (3 - 1.55 ± 0.04 for the high redshift 
sample would result in a over four times higher density param- 
eter Qb(Lya) = 0.044 ■ T'^-''^. Due to the sensitivity of Schayes 
method to yS, the individual treatment without using approxima- 
tions for the column density distribution is to be preferred. 

The unknown gas metallicity does not pose a big error 
source. Though the main metallicit y of the IGM is suspected 
to be about 0.01 solar, according to iCen et al.l (l2005h . most of 
the low column density Lyman a absorbers should have lower 
metallicities (~ lO^^* sol ar). We create d photoionisation models, 
using the Cloudy code dFerla nd et al. 1998) to study the influ- 
ence of metals on the baryon density. The resulting density pa- 
rameter is constant for metallicities lower than 0.01 solar. We 
repeated the calculations with different spectra for the ionising 
background. Using a modified Haardt & Madau spectrum where 
the helium edge is shifted to 3 Ryd, possibly due to the o pac- 
ity of higher He ii Lyman series lines dFechner et al.ll20()6h or a 
background radiation mainly created by AGN as proposed by 
[Mathews & Ferland (1987), resulted in density parameters that 
are within each others Icr error intervals. Misinterpretations in 
the sample of narrow Lyman a lines could add another uncer- 
tainty. This is expected to be small compared to the other men- 
tioned error sources since it is only important for weak absorbers 
which do not contribute much to the total baryon density. 

4.2. BLA detections and WHIM baryon content 

Since the temperature and thus the baryon content of BLAs 
is just determinable when pure thermal broadening can be as- 
sumed, the shape of each broad absorber had to be analysed with 
care. To avoid too many misinterpretations an average signal to 
noise ratio of {S /N) > 15 per pixel was used as a lower limit to 
the data quality. We selected each broad absorber with a Doppler 
parameter of b > 40kms"' by a careful by eye inspection of 
the spectra. Every absorber that showed signs of asymmetries 
in the line profile was excluded from the list. Furthermore no 
saturated absorbers, blended lines or features that are part of 
multi-component systems were used, since the real line width 
and shape is hard to determine in these cases. Since the con- 
tinuum placement can have a large impact on the correct deter- 
mination of the column density and the Doppler parameter, in 
particular for broad and shallow absorbers, just features were in- 
cluded that lie in regions with a clearly defined continuum. 

All thus found absorbers were fitted with the minimising 
algorithm Candalf by R. Baade. The code does a Doppler profile 
fit to the line shape and simultaneously a polynomial fit to the 
continuum. The errors from the line fits are internally added to 
the continuum placement errors. 

It has been empirically found that absorbers with j > 

^^cm /kms" are detectable. Using the minimum signal 

to noise ratio {S /N) = 15 would yield a limit of log j := 

log(j ■ ^JjT^") - 11-1 foi' absorbers that could be detected in 
the whole analysed wavelength range. 

The final sample was divided into two classes. All ab- 
sorbers that still show slight signs of a non-Gaussian shape, 
multi-component substructure, line blends or whose shape is not 
clearly discernable because of noise were marked as tentative. 
At a total redshift interval of Az - 0.723, 38 broad absorbers 
that are possibly originated in the WHIM were found, includ- 
ing 29 tentative cases. Of the 9 good detections, 8 had a sensi- 
tivity coeflicient of log ^ > 11.1 and would thus be detectable 



Table 2. WHIM candidates 



z 


log^H;! 


b [kms-'] 




tent^ 


HE 0515-4414 


1.0624 


13.39 


± 0.08 


51+9 


1.58 


10' 


s 


1.0780 


13.13 


± 0.14 


58 ± 18 


2.07 


10"* 


s 


1.1146 


13.33 


± 0.06 


51 ± 9 


1.59 


10"^ 


s 


1.1276 


13.50 


± 0.08 


101 ± 17 


6.15 


10*^ 


B 


1.1599 


13.44 


± 0.09 


42 ± 8 


1.08 


10' 


SG 


1.1632 


14.02 


± 0.04 


57 ± 5 


1.97 


10' 


s 


1.1884 


13.26 


±0.19 


93 ±42 


5.28 


10' 


B 


1.2665 


13.62 


±0.05 


58 ±5 


2.03 


10' 




1.3206 


13.11 


±0.10 


44± 10 


1.16 


10' 


SG 


1.3395 


13.85 


±0.03 


53 ±4 


1.70 


10' 


BG 


1.5377 


13.41 


±0.05 


41 ±5 


1.04 


10' 




1.5390 


13.83 


±0.02 


54 ±3 


1.74 


10' 




1.5454 


13.01 


±0.26 


108 ± 46 


7.13 


10' 


S 


1.6331 


13.13 


±0.10 


145 ± 22 


1.28 


10*^ 


G 


1.6501 


13.29 


±0.17 


239 ± 53 


3.45 


10« 


GM 


1.6547 


13.19 


±0.19 


177 ± 37 


1.89 


10* 


G 


1.6830 


12.43 


±0.09 


47 ±9 


1.34 


10' 




HE 0141-3932 


1.5945 


12.75 


±0.05 


42 ±5 


1.06 


10' 




1.6000 


12.77 


±0.06 


53 ±6 


1.70 


10' 


S 


1.6168 


12.61 


±0.14 


64± 15 


2.51 


10' 


S 


1.6195 


12.76 


±0.05 


48 ±6 


1.40 


10' 


SB 


1.6285 


12.80 


±0.07 


73 ± 11 


3.20 


10' 


BG 


1.7176 


12.87 


±0.03 


54 ±4 


1.75 


10' 


G 


1.7416 


13.03 


±0.03 


67 ±4 


2.76 


10' 


B 


HE 0429-4901 


1.8607 


12.99 


±0.15 


59± 17 


2.12 


10' 


B 


1.8642 


13.10 


±0.04 


42 ±4 


1.06 


10' 




HE 2225-2258 


1.6926 


13.12 


±0.11 


77 ± 13 


3.56 


10' 


B 


1.7608 


12.54 


±0.13 


50± 13 


1.51 


10' 


B 


1.8031 


12.35 


±0.15 


53± 17 


1.69 


10' 


G 



PG 1634+706 



0.8876 


12.97 


±0.09 


48 ± 


11 


1.38 


10' 


M 


0.9249 


13.02 


±0.07 


53 ± 


8 


1.69 


10' 


SG 


0.9852 


12.75 


±0.19 


67 ± 


26 


2.73 


10' 


GM 


1.1618 


13.86 


±0.02 


43 ± 


1 


1.14 


10' 




1.2038 


13.00 


±0.06 


50 ± 


7 


1.51 


10' 


SG 


1.2255 


13.69 


±0.02 


46 ± 


3 


1.27 


10' 




1.2269 


13.81 


±0.02 


80 ± 


4 


3.83 


10' 




1.2838 


13.39 


±0.06 


41 ± 


5 


1.04 


10' 


B 


1.2846 


13.31 


±0.07 


60 ± 


11 


2.17 


10' 


B 



« logiVm :=logO) 

* The temperatures were determined under the assumption of pure 
thermal broadening. 

^ The tentative detections are marked: (S) 5/A' too low, (B) line 
blend, (G) non-Gaussian shape, (M) multi-component substructure 



throughout the analysed wavelength range. With a total number 
of 9 (38) absorbers, the number of detections per unit redshift 
is ^ = 13.8 (51.2). The bracketed values represent the whole 
sample, including the tentative detections. Table l4!2l shows the 
column density A^hi, the Doppler parameter b and the tempera- 
ture T of all selected broad absorbers. The tentative detections 
are marked in the last column. The most reliable WHIM detec- 
tions are shown in Fig. 1. 
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500 



velocity offset [km s ^ ] 

Fig. 1. Plots of the most reliable WHIM detections sorted by red- 
shift. The Doppler parameter of each absorber is given in km s ' . 
HE 05 15-4414: d, e, f, h. HE 0141-3932: g. HE 0429-4901 : i. PG 
1634-H706: a, b, c. 



The comoving path length was determined using Eq. |2l Just 
parts of the spectrum with a clearly defined continuum, that 
are thus applicable for WHIM detection, were included. The 
resulting total blocking corrected path length has a value of 
AX - 1 .964. Under the assumption of thermal broadening, the 
ionisation fraction was calculated with Eq. |2]for each absorber 
The baryon density was then estimated with Eq.[T] 

The resulting density parameter is Qb(BLA) - (0.91 ±0.05)- 
10"^ for the good sample and Qb(BLA)tent = (8.4 ± 1) ■ 10"^ for 
the whole sample, including the tentative detections. In the lower 
redshift range 0.9 < z < 1.5 a density parameter of Qb(BLA) - 
(1.50±0.07)- 10-^ (Qb(BLA)tent = (4.8 ± 0.6) ■ lO"^) was found. 
The higher redshift interval 1.5 < z < 1.9 yielded Qb(BLA) = 
(0.33+0.02)- 10-3 (Ob(BLA)tent = (12 ± 1) ■ lO^^). Applying the 
strict detection limit of log j > 11.1 resulted for the low redshift 
range in a baryon density of f2b(BLA)tent = (3.4 + 0.4) ■ 10"^ and 
for the high redshift range in Ob(BLA)tent = (0.68 ± 0.07) ■ 10 ^ 
using the whole sample. 

To consider effects of photoionisation the modified ioni- 
sation fraction (Eq. |8]l was used. Thus a baryon density of 
Qb(BLA) = (2.2 ± 0.1) ■ 10-3 (Qb(BLA)te„, = (14 ± 2) ■ lO^^) 
was found. 



5. Discussion 

The density parameter of the total baryonic matter, determined 
from the nucleosynthesis model and measurements of primor- 
dial element abundances, is expected to be Obm ~ 0.04, us- 
ing a Hubble constant of Hq = 72kms-'. A baryon density of 
Qb(LyQ') = 0.010 + 0.001 ■ T^^^ for the diffuse photoionised 



IGM thus represents a fraction of at least ~ 25% to the total 
baryons for temperatures T > 10"* K. The WHIM was found 
to contribute with at least ~ 2.3% (complete sample: 22%) to 
the total baryonic matter Using the modified ionisation fraction 
(Eq. [8]) yielded a contribution of leastwise ~ 5.5% to the total 
baryons. 

While the complete sample should overestimate the baryon 
density, the good sample is expected to underestimate the baryon 
content significantly and can thus be seen as a lower limit. Many 
absorbers that might trace the WHIM were excluded, since non- 
thermal broadening may play a role. Very broad absorbers with 
Doppler parameters b > 100 km s ', corresponding to tempera- 
tures of r > 6 ■ 10^ K, were all marked as tentative since the line 
shape is not discernable with the current data quality available. 
Smallest fluctuations in the quasar continuum make a determi- 
nation of the line shape impossible. Since these very broad and 
thus hot absorbers have a high ionisation fraction, they might 
inhabit large amounts of baryons that remain undetected. Even 
with higher resolution data with a better S /N it could turn out to 
be impossible to trace gas at temperatures T > 10* K with this 
method since the quasar continuum is not completely flat and its 
formation is largely unknown. Effects of non-thermal broaden- 
ing and continuum fluctuations would not be distinguishable. 

Hydrodynamical simulations are able to reproduce the fea- 
tures of the Lyman a forest quite w ell and are thus seen as 
a reliable tool for studying the IGM. I Pave et all (1200 lb found 
that at z a; 1.5 about 70% - 80% of the baryons reside within 
the diffuse photoionised IGM, which is higher than our result. 
The reason can be found in the uncertain intensity of the ionis- 
ing background radiation, the unknown temperature and a pos- 
sibly incomplete line sample. The WHIM is supposed to contain 
10% - 20% of the total baryons. As suspected our result lies 
beneath that value. 

The result indicates that the IGM is in a stage of evolution 
at z K 1.5. While for temperatures T > lO'^K at 1.5 < z < 1.9 
still over 30% of the baryons reside in the diffuse photoionised 
IGM, flie value shrinks to 20% for 0.9 < z < 1.5. The obvi- 
ous increase in the baryon density of the warm-hot phase of the 
IGM from ~ 0.8% for the high redshift sample to ~ 4% for the 
low redshift sample shows that indeed a considerable fraction 
of the baryons has been shock heated to higher temperatures at 
the analysed redshift range. However, the result has to be treated 
with care since the two redshift intervalls were traced with dif- 
ferent instruments, the atmospheric cut-off at /Icut-off ~ 3000 A 
corresponds to a redshift of Zcut-off ~ 1 -5 for hydrogen absorp- 
tion. Thus the low redshift data were taken with STIS, while the 
higher redshift data were obtained with UVES. However, since 
there are more regions with a high S /N in the UVES data then in 
the STIS data, more very broad BLAs that contain high amounts 
of baryons would be detectable in the high redshift range. Thus 
the general trend is not suspected to be created by the change of 
instruments. Though a lower S /N could result in a higher num- 
ber of falsly as broad features interpreted line blends it also in- 
creases the number of lines that are neglected because the shape 
could not be determined satisfactorily. Additionally a systematic 
effect in the change of instruments would have the same influ- 
ence on the narrow as on the broad lines, which is not supported 
by the observations. We found WHIM absorbers up to a redshift 
of z = 1.863. Using the strict detection limit of log j > 11.1, the 
complete sample also shows a considerable increase of baryons 
with lower redshift. The reason for the opposite behaviour of 
the complete sample for the general detection limit lies again in 
the different quality of the UVES and STIS data. Very broad ab- 
sorbers, which contain the highest amounts of baryons, are all 



6 



N. Prause et al.: The baryon density at z = 0.9 - 1.9 



marked as tentative and are just detected in the highest quahty 'PG 1634+706' on paged 
regions of the UVES data. They would not be detectable in the 
STIS data, and thus at lower redshifts. 

It is evident that the WHIM contributes to an important de- 
gree to the total baryon density, not only in the local universe 
but also at higher redshifts 1 < z < 2. A study of the detailed 
evolution of the WHIM, to solve the missing baryon problem 
once and for all, is of great interest. The future awaits several 
X-ray missions that should be able to detect diffuse emission 
from WHIM sources. Until now the search for broad Lyman a 
absorption seems to be the most promising approach. High res- 
olution spectra with a good S /N would be required to make a 
reliable estimation of the baryon content of the WHIM at least 
for the temperature range T ~ 10^ - 10^ K. The Cosmic Origins 
Spectrograph (COS), a new instrument to be installed on the 
Hubble Space Telescope during the 2008 servicing mission, will 
extend the cuiTent STIS data for the low redshift range z < 1 .5. 
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